General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




/ 


* H 


(NASA-CR-158153) THE VISCOELASTIC BEHAVIOR 
OF THE PRINCIPAL COMPLIANCE MATRIX OF A 
UNIDIRECTIONAL GR A PHIT E/EPC X Y COMPOSITE 
Interim Report (Virginia Polytechnic Inst, 
and State Uni*.) 28 F HC A03/MF A01 


N79-16917 


G3/24 


Onclas 

13963 


COLLEGE 






College of Fngineering 

Virginia Polytechnic Institute and State University 
Blacksburg, VA 24061 


VPI-E-79-9 February 1979 

THE VISCOELASTIC BEHAVIOR OF THE 
PRINCIPAL COMPLIANCE MATRIX OF A 
UNIDIRECTIONAL GRAPHITE/EPOXY COMPOSITE 

by 

D. H. Morris^ \ Y. T. Yeow^, and H. F. Brinson^ 


^Department of Engineering Science and Mechanics 
Virginia Polytechnic Institute and State University 
Blacksburg, VA 24061 

^ ^ Allied Chemical Corp. 

Morristown, NJ 07960 


Prepared for: 

National Aeronautics and Space Administration 
Grant No. NASA-NSG 2038 



Materials and Physical Sciences Branch 
Ames Research Center 
Moffett Field, CA 94035 


Approved for Public Release; distribution unlimited. 


ABSTRACT 


A testing program was conducted to determine the time- temperature 
response of the principal compliances of a unidirectional graphite/ 
epoxy composite. It Is shown that two components of the compliance 
matrix are time and temperature Independent. In addition, the compliance 
matrix is found to be symmetric for the viscoelastic composite. 

The time- temperature superposition principle Is used to determine 
shift factors. It is shown that shift factors are independent of fiber 
orientation, for fiber angles that vary from 10° to 90° with respect 
to the load direction. 
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INTRODUCTION 


Composite materials are finding Increasing use In the aerospace 
and automotive Industries. These materials may exhibit time-dependent 
mechanical behavior, depending upon the state of stress, temperature, 
and relative humidity [1]. Thus, It Is necessary to measure the time- 
dependent effects If the materials are to be used for structural 
applications. That such time dependent effects are Important has been 
previously demonstrated [2]. The result Is shown in Fig. 1, which 
indicates that a delayed viscoelastic fracture process was observed 
for a graphite/epoxy [i45] 4$ tensile specimen containing a circular 
hole. That Is, the laminate eventually fractured even though the 
applied load was relaxing in a fixed grip situation. The reason this 
is possible is directly related to the construction of the laminated 
composite being tested. That Is, cracks can open and propagate in one 
ply leading to creep within that ply. This In turn leads to a transfer 
of load to adjacent plies. Thus, while the overall laminate may be 
relaxing, individual fractured plies may exhibit a creep response. 
Obviously, should the same phenomena occur in a prototype structure, 
premature failures would occur. 

It would be desirable to measure the time-dependent effects with 
short-term laboratory tests rather than perform long-term prototype 
studies. It is then clear that there is a need for accelerated 
characterization techniques for composite laminates similar to those 
used for other structural materials. 
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For metals and polymers a variety of accelerated characteriza- 
tion techniques are available such as linear elastic stress analysis, 
empirical extrapolative equations such as the Larson-MIller parameter 
method. Minor's rule and frequency Independence, the time -temperature 
superposition principle, etc. The approach taken In the work reported 
herein Is based upon the time- temperature superposition principle 
developed for polymeric materials, and the widely used lamination 
theory developed for composite materials. 

The procedures for accelerated characterization and lifetime pre- 
diction are outlined In Figs. 2 and 3. Figure 2 Illustrates the plan 
to determine the modulus master curve for unidirectional laminates 
from short-term (15 min.) tensile creep tests from room temperature to 
about 30°C above the glass transition temperature (180°C). Testing 
should occur for various fiber angles from 0° to 90° with respect to 
the load direction. From this series of tests shift functions versus 
temperature and fiber angle are to be determined, as shown In Fig. 2b. 
Linear viscoelasticity will be assumed, and lamina tensile and shear 
strength master curves will be obtained by assuming that a strength 
master curve has the same character as a modulus master curve. 

Probably the most important aspect of the accelerated characteri- 
zation plan Is the generation of a shift function relation, such as 
shown in Fig. 2d. A more detailed discussion of the shift function 
relation will be given later. 

The authors envision the final accelerated design process as 
illustrated by the flow chart shown in Fig. 3. The success of the 
plan shown In Fig. 3 is dependent upon several assumptions. Two critical 
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assumptions occur In boxes B and F of the flow chart of Fig. 3. These 
assumptions Involve symmetry of the principal compliance matrix used 
In modulus transformation equations (box B), and equality of shift 
function with fiber angle (box F). The major thrust of the remainder 
of this paper Is directed toward a critical analysis of the two 
assumptions. Once the assumptions are proved valid, then the remaining 
portion of the plan shown In Fig. 3 should follow and be reasonably 
accura te. 

In order to see how the assumptions fit Into the flow chart of 
Fig. 3, It Is first necessary to establish the constitutive theory for 
a lamina. Assuming a state of plane stress In an orthotropic material, 
the constitutive equation may be written [3], 
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where the x-y coordinate system Is as shown in Fig. 4. If the coordinate 
system is aligned with the lamina principal axes 1-2 (1 denotes the 
fiber direction and 2 denotes the normal to the fiber direction), Eq. 

( 1 ) becomes 
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In Eqs. (1) and (2) the S^- and . are compliance matrices and are 
related through tensor transformation equations [3]. For example, 
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Sj-| * ♦ 2m 2 n^S 12 ♦ n 4 S 22 + m2,l2s 66 (3) 

where m ■ cos 0 and n ■ sin 0 (see Fig. 4). Similar relations exist 
between other components of the compliance matrices [3]. All primed 
quantities are only a function of the four principal compliances of 
Eq. (2) and the angle 0. 

Thus, one need only determine the principal compliances from Eq. 
(2) to completely define the state of stress or strain for a continuous 
fiber reinforced composite under a state of plane stress. If the 
material Is viscoelastic, the components of the principal compliance 
matrix will be a function of time, temperature, relative humidity, and 
stress level [1,4,5]. The present study will consider stress levels 
such that material response Is linear viscoelastic. In addition, 
relative humidity will not be an experimental variable. 

When the composite material response Is linear viscoelastic, 

Eq. (3) may be written [5] 

Sj^t) - m 4 S n (t) + 2m 2 n 2 S 12 (t) + n 4 S 22 (t) ♦ m 2 n 2 S 66 (t) (4) 

where t denotes time. Equation (4), together with composite lamination 
theory, may be used to predict the viscoelastic behavior of a general 
laminate. 

In brief, the plan consists of constructing master compliance 
curves for unidirectional laminates from short-term creep tests at 
various temperatures. This requires the use of the time- temperature 
superposition principle. From the short-term tests the shift function 
versus temperature relationship may also be found. Using the master 
compliance curves, the shift functions, and equations similar to 
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Eq. (4), compliance master curves for any angle e (Fig. 4) could be 
constructed. These master curves could then be used as Input to an 
Incremental computational procedure using standard stress analysis 
lamination theory to predict long-term material behavior. 

As previously mentioned the above procedure Is dependent upon the 
use of transformation equations, such as Eq. (4), and the shift 
function- temperature relationship. Two comments are In order regarding 
the previous statement. First, In writing Eq. (4/ It Is assumed that 
the compliance matrix In Eq. (2) Is symmetric, that Is, 

Si 2 (t) ■ S 21 (t). And seconi, for Eq. (4) to be of practical use, the 
shift functions for the components of the principal compliance 
matrix, Eq. (2), should be equal. 

Based on these corments, the objective of this paper is twofold; 
(1) check the assumption that 2 ( t ) * S^(t), and (2) show whether 
the shift function-temperature relationship is the same for various 
angles 9, that Is, show whether the shift function Is Independent of 
fiber orientation. 


EXPERIMENTAL CONSIDERATIONS 

The particular composite material studied In this investigation 
was manufactured from T300/934 graphite/epoxy pre-preg tapes. All 
test specimens were cut from a single large panel of the material. In 
addition, all specimens were unidirectional laminates. 

Tensile specimens were instrumented with three-element rectangular 
strain gage rosettes. Two rosettes were bonded to each specimen, one 
on each side, and the gage outputs were averaged to eliminate any 
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out-of-plane bending. The strain gage rosettes were temperature 
compensated using specimens whose fiber orientation was Identical to 
the stressed specimens. 

load was applied to test specimens through special grips, similar 
to those used by Chamls and Sinclair [6]. They show that these grips, 
together with the specimen length used, tended to minimize the In-plane 
bending discussed by Pagano and Hal pin [7]. The grips used by Chamls 
and Sinclair were modified by the addition of a pin that extended 
through the midpoint of the clamped specimen section. This pin helped 
to reduce slippage between the metal clamps, bonded end tabs, and test 
specimen. Load levels were such that material response was always 
linear viscoelastic. Further discussion of the experimental procedures 
may be found elsewhere [8]. 

ANALYTICAL CONSIDERATIONS 

As previously stated, one can completely characterize lamina 
viscoelastic behavior by determining the principal compliances S^(t), 
S 22 (*)» and Sgg(t) that appear In Eq. (4). Thus, In order to 

answer the questions raised In the objective statement. It Is appropriate 
to discuss the methods whereby the components of the principal 
compliance matrix were determined. 

First, let 0*0 (Fig. 4) and write the principal compliance 
S n (t) as 

«,(t) 

(5) 

0 

where e^(t) denotes that the axial strain is a function of time, and 
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<3 q Is the constant applied stress. The axial component of strain (as 
well as those that follow) was found using the strain gage rosettes. 

A c.*eep test of the same specimen simultaneously yielded the value of 
S 2 (t) from the expression 

e 2 (t) 

S 12 (t) ' (6) 


where e 2 ( t) Is the time dependent transverse strain. 

In order to check the assumption of symmetry of the principal 
compliance matrix, It Is necessary to experimentally determine S 21 ( t ) 
from 

e,(t) 

S 2 ,(t) . 4 - (7) 

0 


where e ■ 90° in Fig. 4. Note that ( t ) does not have the same meaning 
In Eqs . (5) and (7) . 

The two remaining principal compliances are found as follows. 
Letting 0 *90°, the compliance S 22 U) may be written as 

«,(t) 

S 22 (t) . 4- (8) 


where § 2(0 Is the time dependent axial strain, 
was found using 


* 56 ^ 


Y 12 (t) 
T 1 2 


The compliance S 66 (t) 
(9) 


where 0 a 10°. Further discussion on the applicability of Eq. (9) for 
calculating Sg 6 (t) may be found in Ref. [8]. In Eq. (9), 1* a 

constant intralaminar shear stress and is found via the stress 
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transformation equation 

t 12 * \ °o s ^ n 20 (10) 

The time dependent shear strain, Is determined from the strain 

transformation equation 

y 12 * ' * e x * e y* s1n 26 + Y xy C0S 20 O 1 ) 

where 

y * 2e„ c - e - e (12) 

'xy 45 x y 

In Eqs. (11) and (12), e x> and e 45 are the axial, transverse and 
45° strains, respectively, obtained from the three-element rectangular 
strain gage rosettes. 

Equations (5) - (12) may be used to evaluate the components of the 
principal compliance matrix. In particular, Eqs. (6) and (7) may be 
used to determine the symmetry of the compliance matrix. Using Eqs. 

(5), (6), (8), d-.d (9) and short-term experimental results, one can 
construct master compliance curves and find shift functions, and thus 
determine whether the shift functions are independent of fiber orienta- 
tion, e. 


RESULTS AND DISCUSSION 
Symmetry of Compliance Matrix 

Table 1 shows the experimental results for S-^t) and S 21 ( t ) , 
for t * 1 min. Similar results were found for other values of time up 
to 15 min., which was the extent of the short-term tests. Several 
results are evident from inspection of the data in Table 1. Both S^ 
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and S 21 are essentially Independent of temperature (within about 10%). 

The maximum temperature was 210°C. Data at t * 15 min. indicated that 
and are within about 10% of the values of Table 1, which shows 
that these compliance terms are also essentially Independent of time. 

The average values of S 12 and $ 2 ^ In Table 1 are within 9% of each 
other. The differences cited above are probably due to experimental 
error due to the small transverse strains measured In a e ■ 90° test, 
and due to the fact that only one test was conducted at each temperature 
level. It is felt that the scatter in data represents usual variations 
in test results encountered in composite materials. Henceforth, the 
compliance matrix will be considered symmetric. In addition, it will 
be assumed that independent of time and temperature. Experimental 

results also indicate that is independent of time and temperature. 
Similar results have been found for a unidirectional glass fiDer-epoxy 
composite material [1], 

Equality of Shift Functions 

Before proceeding with a discussion of the equality of shift 
functions, it is appropriate to briefly discuss the methods whereby the 
results were obtained. Further details may be found in [9], 

Figure 5 illustrates the time-temperature behavior of the reduced 
reciprocal of S g6 (e ■ 10°), for short-term (15 min.) creep. The 
reduced compliance was calculated from 

6p(f) T 

s 66 (t) ’-f— r < 13 > 

0 0 

where T represents the absolute test temperature and T Q was the absolute 
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reference temperature (taken as the glass transition temperature of 
453°K, [8]). The ordinate of Fig. 5 was plotted In terms of reduced 
reciprocal of compliance In order to study the applicability of the 
time-temperature superposition principle to composite materials [8,9]. 
Figure 5 also shows a portion of the master curve at the reference 
temperature of 453°K, or 180°C. A complete master curve was obtained 
by graphically shifting the short-term curves along the log time axis 
until one continuous curve was obtained. The shift factor was the 
amount of horizontal shifting necessary to superpose the various 
constant temperature curves. 

Figure 6 shows shift factors, ay, for fiber orientations that 
range from 10° to 90°. For fiber orientations of 15°, 30°, 45°, 60°, 
and 75°, compliances were calculated using an equation similar to Eq. 
(9), and Eqs. (10), (11), and (12). For all practical purposes the 
shift factors in Fig. 6 are equal, and thus are Independent of fiber 
orientation. A similar result has been found by Moehlenpah, Ishai, and 
DiBenedetto [10] for the tensile yield stress shift factors for glass 
fiber reinforced epoxies. Note that and are Independent of 
time and temperature and do not require the use of shift functions. 

As seen in Fig. 6, there is scatter in experimental data. How- 
ever, above the glass transition temperature of 180°C data scatter is 
slight. Some of the reasons for scatter, or lack thereof, may be seen 
by examination of Fig. 5. Above 180°C the slopes of the compliance- 
time curves are greater than for temperatures less than 180°C. Thus, 
graphical shifting was more accurate for temoeratures above the glass 
transition temperature. In addition, as previously stated, only one 
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test was conducted at each temperature level. It Is felt that the 
data scatter Is typical of that found when testing composite materials. 
For more accurate results several tests at a particular temperature 
should be conducted. 


CONCLUSIONS 

It has been shown that the principal compliance matrix for an 
orthotropic body under plane stress Is symmetric. Two components of 
the matrix, S^ and S^» are time and temperature independent. Similar 
results have been found for a glass-epoxy unidirectional composite [1]. 
In addition, shift factors are independent of fiber orientation for 
fibers oriented between 10° and 90° to the load direction. 

Both of the above results are important when using the tensor 
transformation equations to predict the time-dependent compliant be- 
havior of off-axis unidirectional laminates, and hence general laminate 
behavior as shown in the plan of Fig. 3. 

Obviously, if reliable values of principal compliances, master 
curves, and shift factors are to be obtained, then many tests will be 
needed to establish the statistical variation of properties with time 
and temperature. That was not the purpose of the work reported herein. 
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Table 1 . 


Measured Properties and S 21 (1 min.) for 
T300/934 Graphite/Epoxy Composite. 


Temperature Si? x ]0 6 psi" 1 S?i x 10 6 psl* 1 

(°C) (x 10 4 MPa-1) (x 104 MPa-1) 


22 

-0.0149 

-0.0143 


(-0.0216) 

(-0.0207) 

100 

-0.0124 

-0.0148 


(-0.0180) 

(-0.0215) 

180 

-0.0132 

-0.0144 


(-0.0191) 

(-0.0209) 

200 

-0.0140 

-0.0147 


(-0.0203) 

(-0.0213) 

210 

-0.0136 

-0.0157 


(-0.0197) 

(-0.0228) 



15 



125 


100 


75 


50 


25 


Fig. 1. Time-Dependent Fracture of [+45L T300/934 Graphite/Epoxy 

Plate With a Central Circular Hole 
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a) MEASURED TENSILE MODULUS 
MASTER CURVES VS. FIBER 
ANGLE FOR UNIDIRECTIONAL 
LAMINATE. 



b) MEASURED SHIFT FUNCTION VS. 
FIBER ANGLE FOR UNIDIRECTIONAL 
LAMINATE. 



c) TENSILE STRENGTH MASTER 
CURVE ASSUMED TO HAVE 
SAME CHARACTER AND SHIFT 
FUNCTION AS MODULUS 
MASTER CURVE (^MEASURED 
DATA ) . 


d) FORMAL SHIFT FUNCTION VS. 

FIBER ANGLE AND TEMPERATURE 
RELATION ANALOGOUS TO WFL 
EQUATION DEVELOPED FROM b). 


Fig. 2. Proposed Method to Determine Tensile Strength Master Curve 
and Shift Function Equation 
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Fig. 3. Flow Chart of the Proposed Procedures for Laminate Accelerated 
Characterization and Failure Prediction 







l/s xx x 10 (MPa) 


19 


A 23*C, 115*0,175*0, 204*C 
□ 45 W C, I33*C, I80*C, 2IO # C 
O 78 # C, I50*C, 185° C 
V 90 # C, I60*C, I93 # C 
O I0I*C, I70*C, I98*C 



LOG TIME (minutes) 


Fig. 5. Reduced Reciprocal of Compliance and Portion of 180°C Master 
Curve for [10°]g s T300/934 Graphite/Epoxy Laminate 
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Fig. 6. Shift Factors versus Temperature for Creep of Off-Axis Tensile 
Coupons (T300/934 Graphite/Epoxy) 




